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Abstract 

The Kadomtsev-Petviashvili II (KPII) equation admits a large variety of multi-soliton solutions which 
exhibit both elastic as well as inelastic types of interactions. This work investigates a general class of 
multi-solitons which were not previously studied, and which do not in general conserve the number of 
line solitons after interaction. The incoming and outgoing line solitons for these solutions are explicitly 
characterized by analyzing the T-function generating such solutions. A special family of A^-soliton solu- 
tions is also considered in this article. These solutions are characterized by elastic soliton interactions, in 
the sense that amplitude and directions of the individual line solitons as y ^ oo are the same as those of 
the individual line solitons as y ^ — oo. It is shown that the solution space of these elastic A^-soliton solu- 
tions can be classified into {2N — 1 ) ! ! disjoint sectors which are characterized in terms of the amplitudes 
and directions of the line solitons. 



1 Introduction 

The purpose of this work is to present a characterization of a large family of real, non-singular, line-soliton 
solutions of the KPII equation 

{-4Ut + UxXX + 6UUx)x + ^Uyy =0 (1.1) 

where u = u{x,y,t) and subscripts x, y and t denote partial derivatives. The KP equation is perhaps the 
prototypical (2-1-1 )-dimensional integrable evolution equation originally derived ifTOl as a model for small- 
amplitude, long-wavelength, weakly two-dimensional solitary waves in a weakly dispersive medium, and 
arises in many different applications including water waves and plasmas (for a review, see e.g. fTl). The 
various aspects related to the integrability of the KP equation have been studied extensively, and a large 
number of exact solutions have been found. These works are documented in several monographs (see e.g., 
Refs. |[il|7j[l7l and references therein). There are two versions of the KP equation depending on the sign of 
the dispersion namely, KPI (for positive dispersion) and KPII (for negative dispersion). Here we consider 
the KPII equation. 

Among the exact solutions of KPII, perhaps one of the most well known class of real, non-singular 
solutions is the line-soliton solutions. The simplest type is the one-soliton solution which is a traveling wave 
in xj-plane, and is localized along a line. A straightforward generalization of this solution to a multi-soliton 
configuration of line solitons was also found in earlier works Il20l l5l. In the generic case, these A^-soliton 
solutions form a pattern of intersecting straight lines in the xy-plane apart from small spatial shifts arising 
from the pairwise interactions between any two lines. We refer to these solutions as the ordinary A/^-soliton 
solutions which can be parametrized by 2N parameters namely, amplitude and direction (in the xy-plane) of 
the A'^ line solitons. However, it has been shown theoretically and experimentally that it is not possible to 
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obtain ordinary A/^-soIiton solutions for all choices of the soliton parameters (see e.g., Ref. H)- That is, the 
parameter space of ordinary A/^-soliton solutions is not simply the A/^-fold Cartesian product of the parameter 
space of one-soli ton solution of the KPII equation. Subsequent work has revealed that other families of 
line soliton solutions exist in addition to the ordinary A/^-soliton solutions. The simplest of these solutions 
describes the so called "Y-junction" solutions, which describes resonant interaction of two line solitons, and 
has been known since 1977 lITSl [T6l . In this case, the three solitons with wave-numbers and frequencies 
{ka,0)a), a = 1,2,3, satisfy the three-wave resonance conditions: ki +k2 = ks and coi + 0>2 = COj. More 
general resonant solutions have also been obtained in Refs. lfT4l[T8l [T9l. Recently, in Ref. 0, the authors 
found a large family of soliton solutions in which an arbitrary number N- of incoming line solitons interact 
resonantly via intermediate line solitons which produce web like patterns in the xy-plane, and then form an 
arbitrary number A'^_|_ of outgoing line solitons, where ^ N^, in general. We refer to these multi-solitons 
solutions as the (A^-jA/^^) -soliton solutions of KPII. In particular, the case = = N yields yet another 
kind of A^-soUton solutions which however differ significantly from the ordinary A/^-soliton solutions in their 
interaction patterns. The current work is motivated by these recent studies which indicate that the solitonic 
sector of the KPII equation is richer than previously thought as many other families of line-soliton solutions 
exist in addition to the ordinary A/^-soliton solutions. 

In this article we are primarily concerned with the characterization of the incoming and outgoing line 
solitons of a generic multi-soliton configuration of KPII as well as the classification of a particular class 
of multi-soliton solutions called the elastic A/^-soliton solutions. The paper is organized as follows. In 
section [2l we discuss the asymptotic behavior of the T-function underlying the multi-soliton solutions of 
KPII. In particular, we show that as j — > ±oo, the solution u{x,y,t) decays exponentially in the x3'-plane 
except along certain rays which correspond to the incoming and outgoing line solitons. We then characterize 
these asymptotic structures in terms of the parameters of the T-function. In section[3l we study a special class 
of solutions called the elastic A/^-soliton solutions. We show that these solutions can be classified into (2A'^ — 
I)!! inequivalent types determined by the 2A^ soliton parameters comprising of N pairs of amplitudes and 
directions associated with the incoming or outgoing line solitons. Moreover, from a given set of admissible 
soliton parameters one can explicitly determine an equivalence class of elastic A/^-soliton solutions such 
that any two solutions in the equivalence class exhibit similar interaction patterns and have the same set 
of incoming and outgoing line solitons. In this paper, we only present the results and discuss some of 
the important features of the multi-soliton solutions of KPII. The proofs of these results and more detailed 
discussions can be found in Refs. SB!. 

We point out that the elastic A^-soliton solutions of KPII equation was also recently addressed in Ref. fTPl 
where an elegant characterization of these solutions were presented in terms of the Schubert cell decompo- 
sition of the Grassmann variety Gr{N,2N). Here we follow a different approach motivated by the physical 
problem of identifying the distinct types of elastic A/^-soliton solutions with the the corresponding parameter 
space of soliton amplitudes and velocities. Finally, we note that line soliton solutions with novel web like 
spatial structures were also found recently in several other 2 + 1-dimensional integrable equations. Examples 
include Refs. [8, ^ for a coupled KP system, and Ref. |[T2l where similar solutions were found in discrete 
integrable systems such as the two-dimensional Toda lattice and its fully- and ultra-discrete analogues. 

2 Asymptotic line solitons 

In this section we investigate the line solitons of the KPII equation and the asymptotic properties of the 
T-functions generating such solutions. The solution u{x,y,t) of the KPII equation can be obtained from the 
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T-function T{x,y,t) via the relation 



u{x,y,t) =2{\ogT{x,y,t))j,x.i (2.1) 
It is well-known (see, e.g. Refs. ||5l|6l[T3l) that T{x,y,t) can be expressed in the form of a Wronskian 



t(x,3;,0 =Wr(/i,...,/;v) =det 



/ /i /2 ••• In \ 
f\ fi ■■■ In 



(2.2) 



'2 ./W / 

where = fn/dx^ and where the functions form a set of linearly independent solutions of the 

linear system 

fy — fxx 1 ft — fxxx ■ 

A general family of multi-soliton solutions can be constructed in a simple way from Eq. (12.21 ) by choosing 
each function f„{x,y,t) to be a linear combination of real exponentials. That is, 

/n(-^,>',0 = I "', n = l,2,...,N, (2.3) 

m=l 

where, 6^ = k,„x + kf„y + ^J]/ + 0om for ni = 1 , . . . ,M are M phases with real phase parameters ki,...,kM 
and real constants 0oi , . . . , 0om> and where the constant coefficients <3„,„ define the N x M coefficient matrix 
A := (anm)- Note that one can naturally identify each /„ with the row and each phase d,„ with the 
column of the coefficient matrix A, and vice versa. Upon substituting Eq. (12.31 ) into the Wronskian 
of Eq. (12.21 ) and then using the Binet-Cauchy formula to expand the resulting determinant, we obtain the 
following explicit form of the T-function 

T{x,y,t)= A(mi,...,mA,) exp[0(mi,...,mA?)] Yl (^m,-^mj, (2-4) 

l<mi<-<mN<M l<s<r<N 

where A(/«i , . . . , m^) is the N x N minor of A obtained by selecting the columns 1 < mi < • • • < nif^ < M, 
and 6 {mi,..., mj^) := + • • • + ^mjv is a phase combination of N (out of M) distinct phases. The T-function 
given above could in general, vanish at points {x,y,t) € where the solution u{x,y,t) in Eq. (12.11 ) would 
then have singularities. However, the following restrictions on the phase parameters {k„}^^i and on the 
coefficient matrix A are sufficient to guarantee that the resulting solutions u{x,y,t) of the KPII equation are 
nonsingular. 

Condition 2.1 (Positive definiteness) 

(a) The phase parameters are distinct. Hence, without loss of generality, they can be ordered as k\ <k2< 
...<kM. 

(b) The N xM coefficient matrix A satisfies rank(A) = A^, and M > N. 

(c) All non-zero N xN minors of A are positive. 

From Condition 12. H it is clear that the coefficient of each exponential term of the sum in Eq. (12.41 ) is positive 
because the phase parameters ki,... ,kM are well-ordered and all the minors A(mi, . . . ,mN) of A are also 
nonnegative. As a result, T{x,y,t) is a nonvanishing, positive function for all {x,y,t) G R^, and generates a 
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nonsingular solution of the KPII equation via Eq. (12.11) . If M < N, T{x,y,t) = because the set of functions 
{fn}^=i in Eq. (12.21 ) is hnearly dependent; also when M = N, T{x,y,t) in Eq. ( 12.41 ) contains only one expo- 
nential term which leads to the trivial solution u{x,y,t) = in Eq. (12.11 ). Therefore, for nontrivial solutions 
we must have M > N when there are more than one exponential term in the sum of Eq. ( 12.41 ). 

The simplest example is a one-soliton solution obtained by choosing A'^ = 1,M = 2 and f{x,y,t) = 
gQi +e^2, with k\ < kj. This choice yields the following traveling-wave solution 



where x = {x,y) and where the wave vector k := {lx,ly) = (^i —k2,kj — k\) and the frequency (O :=k\ — k\ 
satisfy the nonlinear dispersion relation 



For fixed t, the solution u{x,y,t) decays exponentially in the x3'-plane except along the line Q\ = 62 whose 
normal has a slope c = ly/lx = ^1 +/^2- Such solitary wave solutions of the KPII equation are called line 
solitons. Apart from a trivial constant d\ Q — 62,0 in Eq. ( 12.51 ) corresponding to an overall translation, a line 
soliton of KP is characterized by the phase parameters ki,k2, or by two physical parameters, namely, the 
soliton amplitude a:=k2 — k\ and the soliton direction c :=ki+ ^2- 

When c = (i.e., k\ = —k2), the solution in Eq. ( 12.51 ) becomes ^-independent and reduces to the one- 
soliton solution of the Korteweg-de Vries (KdV) equation. However, due to the dependence on the additional 
spatial variable y, the multi-soliton solution space of the KPII equation is much richer than that of the KdV. 
Indeed we find that Eq. (12.31 ) with the coefficient matrix A satisfying Condition 12. 1 1 leads to a large class of 
multi-soliton configurations which consist of incoming (i.e., as 3^ — > —00) and N+ outgoing (i.e., as 3^ ^ 00) 
line solitons. The amplitudes, directions and the number of incoming line solitons are in general different 
from those of the outgoing line solitons. In order to characterize the incoming and outgoing line solitons 
associated with these solutions, it is necessary to examine the asymptotic behavior of the T-function in the 
xy-plane as |j| ^ 00, and for finite t. Recall from Eq. (12.41 ) that x{x,y,t) is a linear combination of exponential 
phase combinations with positive coefficients. The leading order behavior of the T-function as j ^ ±00 in 
a given asymptotic sector of the x3'-plane is governed by that exponential term which is dominant in that 
region. The solution u{x,y,t) generated by the T-function is exponentially small at all points in the interior 
of any dominant region, and is localized only at the boundaries of the dominant regions, where a balance 
exists between two or more dominant phase combinations in the T-function of Eq. (12.41 ). The asymptotic 
properties of the faw-function and the solution u{x,y,t) can be derived from a systematic analysis of these 
dominant exponential phases. These are summarized below. 

Proposition 2.2 For finite values oft, and for generic values of phase parameters ki,. . . , kM, the incoming 
and outgoing line solitons of the {N-,N^) -soliton solutions of KPII are characterized as follows. 

( i) As y ^ ±00, the dominant phase combinations of the T-function in adjacent regions of the xy-plane 
contain N —\ common phases and differ by only a single phase. The transition between any two such 
dominant phase combinations d{i,m2, ■ ■ ■ ,mt^) and G{j,m2, ■ ■ ■ ,mi^) occurs along the line defined by 
Lij : Bj = 6j, where a single phase Qi in the dominant phase combination is replaced by a phase Gj. 

( ii) Along the single-phase transition line L^j, the asymptotic behavior of the T-function as y ^ ±00 is 
determined by the balance between the two dominant phase combinations d{i,m2, ■ ■ ■ ,m!^) and 
S (7 5 "^2 5 • • • ) w^A?) in Eq. (12.41) . and is given by 



u{x,y,t) = ^{k2 - kif sech^ ^{62 - di) =(£>{k-x + cot) , 



(2.5) 



4colx + 4 + 3ly=0. 



(2.6) 



x,y,t) ^Cic 




4 




Figure 2.1: Dominant phase combinations in the different regions of the x, y-plane (labeled by the indices in parenthe- 
ses) and the asymptotic line solitons (labeled by the indices in square braces) for two different line soliton solutions: 
(a) a Miles resonant solution (Y-junction) with {ki^k2^kj,) = (— |, — j,2) at f =0; (b) an inelastic (2,2)-soliton with 
{k\,...^kA,) — (— |, — j,0, 1) at t — 0. Here and in all of the following figures, the horizontal and vertical axes are 
respectively x and y, and the graphs show contour lines of logM(x,3',f) at a fixed value of t. 



The coefficients Ct and Cj above depend on the phase parameters ki,kj,km2 ,■■ ■ jhuN' on the N xN 
minors A{i,m2, ■ ■ ■ ,mN) and A{j,m2, ■ ■ ■ ,m!^) of the coefficient matrix A. The asymptotic behavior of 
the solution in a neighborhood of a single-phase transition is then obtained from Eq. (12.11) as 

u{x,y,t) ~ \{ki-kjf^&ch^ [\{Bi-Bj)\ , (2.7) 

which is a traveling wave satisfying the dispersion relation in Eq. (12.61 ). Equation (12.71) has the same 
form as the one-soliton solution in Eq. i2.5\) , and thus it defines an asymptotic ( incoming or outgoing) 
line soliton associated with the single-phase transition i j. For each asymptotic line soliton, the 
soliton amplitude is given by atj = \ki — kj\, and the soliton direction is given by cij = ki + kj, which 
the direction (slope of the normal vector) of the transition line Ljj. 

We can label each asymptotic line soliton associated with the single-phase transition / — > j, by the index pair 
[ij] which uniquely identifies the phase parameters kj and kj in the ordered set {^i, . . . j^m}- 

The simplest instance of a transition of dominant phase combinations arises for the one-soliton so- 
lution (12.51) . which is localized in the xy-plane along the line di = 62 defining the boundary of the two 
half-planes where each of the two phases Bi and B2 dominates. In the general case, the dominant regions 
are more complicated although the solution u{x,y,t) is still localized along the boundaries of these regions. 
For example. Fig. 12. H a) illustrates a (2, l)-soliton solution mentioned in the introduction as a Y-junction 
lITSl describing the resonant interaction of two line solitons. This solution corresponds toN = l,M = 3, and 
is generated by the T-function z{x,y,t) = e^' -\-e^- + e^^. In this case, the xy-plane is partitioned into three 
dominant regions corresponding to each of the dominant phases di, 62 and 63. Once again, the solution 
u{x,y,t) is exponentially small in the interior of each dominant regions, and is localized along the phase 
transition boundaries: di = 62, di = 63 and 62 = 63. Each of the asymptotic line solitons labeled by the in- 
dex pairs [1,2], [2,3] and [1,3] are given by Eq. (12.71) and satisfy the one-soliton dispersion relation Eq. (12.61 ). 
While some of the dominant regions have infinite extensions in the xy-plane, others can be bounded, as in 
the case of the (2,2)-soliton shown in Fig. I2.1f b). This solution is generated by the T-function in Eq. (12.21 ) 
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with /i = e^' - e^^ and /2 = e"^ + e^^ + e""^, i.e., N = 2,M = 4. In addition to the unbounded dominant 
regions corresponding to the phase combinations 0(1,2), 0(1,4), 0(3,4) and 0(1,3), in this case there is 
also a bounded region in the Aj-plane where 0(2,4) is the dominant phase combination. The boundaries of 
this region is formed by the incoming asymptotic line solitons [1,4] and [2,3], together with the intermediate 
line soliton [1,2]. Note that the outgoing line solitons [1,3] and [2,4] interact resonantly via two Y-junctions, 
while the incoming soliton pair interact non-resonantly. Note also that the intermediate segment is a line 
soliton in its own right, since the solution is locally given by Eq. (I2.7I ). with [i,j] = [2,3]. 

According to Proposition [Z21 an asymptotic line soliton corresponds to a dominant balance between two 
phase combinations in the T-function. But one still needs to identify which particular phase combinations 
are indeed dominant in a given T-function as \y\ — > oo. This requires a closer examination at the structure of 
the N X M coefficient matrix A associated with the T-function. Note that elementary row operations on A 
given by A —5- A' = GA where G € GL(A'^,R), amounts to an overall rescaling of the T-function in Eq. (12.41 ). 
i.e., ■T{x,y,t) —>■ ■T'{x,y,t) = det{G)x{x,y,t). Since such rescaling leaves the solution u{x,y,t) in Eq. (12.11 ) 
invariant, it is possible to choose the coefficient matrix A in reduced row-echelon form (RREF) making use 
of Gaussian elimination. Recall that, for anN xM matrix in RREF, the leftmost nonvanishing entry in each 
nonzero row is called a pivot, which is normalized to 1 , so that the pivot columns are the elements of the 
canonical basis of M.^. Throughout the rest of this work we will consider the coefficient matrix A to be in 
RREF, and to satisfy Condition 12. H and the following additional conditions: 

Condition 2.3 (Irreducibility) Each column of A contains at least one nonzero element, and each row of A 
contains at least one nonzero element in addition to the pivot. 

Then a detailed analysis of the coefficient matrix A satisfying Conditions 12.11 and 12.31 leads to an explicit 
identification of those / j single-phase transitions which actually occurs as |y| — > °°, for any given T- 
function of Eq. (12.41 ). As a result, each asymptotic line-soliton [i,j] is also explicitly determined by the 
coefficient matrix A. as given below. Specifically, we have the following results. 

Proposition 2.4 Each index pair [i,j] labeling an asymptotic line soliton of a (N-,N+)-soliton solution are 
uniquely identified with a pair of columns of the associated coefficient matrix A, as prescribed below. 

( i) An asymptotic line soliton as y is identified by a unique index pair [e„ , jn] with e„ < j,, and where 
{^n}^=i label the pivot columns of A. Similarly, an asymptotic line soliton as y ^ —°° is identified 
with a unique index pair [in,gn] with /„ < g„ and where label the non-pivot columns of A. 
Thus, the {N-,N^)-line soliton solution of KPII generated from the T-function in Eq. (12.41 ) has exactly 

= N asymptotic line solitons as y ^ °° and N- = M —N asymptotic line solitons as y ^ —oo. 

(ii) The necessary and sufficient conditions for an index pair [i,j] to identify an asymptotic line soliton 
is determined by considering the ranks of two sub-matrices Xjj and Yij of A. They can be denoted by 
their column indices as follows: 

= [1,2, 1,7 + 1, .. . ,M] Yij = [/+ 1, . . . j- 1] . 

That is, Xij consists of all consecutive columns to the left of the column and all consecutive columns 
to the right of the column of A, while Yij consists of all consecutive columns in between the and 
column of A. The rank conditions are then stated as follows. 

(a) [i,j] identifies an asymptotic line soliton as y ^ oo if and only if rank(X,y) := r <N — I and 
mnk{Xij\i) = rank(X,y | j) = rank(X,y|/,7) = r + 1. 
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(b) [i,j] identifies an asymptotic line soliton as y ^ —°° if and only if rank(y,y) := s <N —\ and 
rank(y,y|/) =rank(F;y|j) = rank(y,-y|/, j) =s+ 1. 

Above, {Z\m,n) denotes the sub-matrix Z of A augmented by the columns m and n of A. 

Note that for the asymptotic Une soliton [e„,i„] as y ^ oo in Proposition I2.4f i). e„ is a pivot index but 
the index ]„ can be either a pivot or a non-pivot index. Similarly, for the asymptotic line solitons [4,^,,] 
as 3^ ^ — oo, the index g„ is an non-pivot index, while /„ can be either a pivot or a non-pivot index. The 
necessary and sufficient rank conditions in Proposition I2.4f ii) provides a constructive method to identify the 
asymptotic line solitons as y ^ ±00 from a given coefficient matrix A in RREF. We illustrate these statements 
by the examples below. 

Example 1 : Consider the T-function in Eq. (12.41 ) with N = 2 and M = 5 generated by the coefficient matrix 

The pivot columns of A are labeled by the indices {e\,e2} = {1,3}, and the non-pivot columns by the indices 
{g\,g2,g3} = {2,4,5}. Thus, from Proposition I2.4f i) we know that there will be N+ = N = 2 asymptotic 
line solitons as y — > 00, identified by the index pairs [l,7i] and [3,72] for some ji > 1 and 72 > 3, and that 
there will be = M — N = 3 asymptotic line solitons as j — > —00, identified by the index pairs [i\ , 2] , [i2 , 4] 
and [/3,5], for some ii < 2, i2 < 4 and 23 < 5. We first determine the asymptotic line solitons as 3^ ^ 00 
using the rank conditions prescribed in Proposition I2.4f ii). For the first pivot column, ei = 1, we start with 
j = 2 and consider the sub-matrix Xi2 = {^ i i) ■ Since rank(Xi2) = 2 which is greater than — 1 = 1, we 
conclude that the pair [1 , 2] cannot identify an asymptotic line soliton as j ^ 00. Incrementing 7 to 7 = 3, 4, 5 
and checking the rank of each sub-matrix Xij we find that the rank conditions in Proposition I2.4f ii) are 
satisfied when 7 = 4: Xu = {^i), so rank(Xi4) = 1 and rank(Xi4|l) = rank(Xi4|4) = rank(Xi4|l,4) = 2. 
Thus, the first asymptotic line soliton as j — 00 is identified by the index pair [1,4]. For the second pivot, 
^2 = 3, proceeding in a similar manner we find that 7 = 4 does not satisfy the rank conditions (since X34 
has rank 2) but 7 = 5 does: X35 = ("j* ^j^) , which yields rank(X35) = 1 and rank(X35|3) = rank(X35|5) = 
rank(X35|3,5) = 2. Therefore, the asymptotic line solitons as j — > 00 are given by the index pairs [1,4] and 
[3,5]. 

Next we consider the asymptotics for y — > —00. Starting with the non-pivot column g\ = 2, the only 
column to its left is / = 1. Then, we have Y\2 = 0, and rank(Fi2|l) = rank(Fi2|2) = rank(yi2|l,2) = 1. 
Consequently, the pair [1,2] identifies an asymptotic line soliton as y — > —00. For g2 = 4 we consider 
/= 1,2,3 and find that the rank conditions are satisfied only for / = 2. In this case, Y24 = (j), sorank(F24) = 
I = N — I and rank(F24|2) = rank(}24|4) = rank(y24|2,4) = 2. Hence [2,4] is the unique asymptotic line 
soliton as J — > —00 associated to the non-pivot column g2 = 4. In a similar way, we can uniquely identify 
the last asymptotic line soliton as y — > —00 as given by the indices [3,5]. To summarize, there are = 2 
outgoing line solitons given by the index pairs [1,4] and [3,5], and there are = 3 incoming line solitons 
given by the index pairs [1,2], [2,4] and [3,5]. A snapshot of the solution at f = —32 is shown in Fig. 12.21 

Example 2: Next, consider the T-function with N = 3 and M = 6 generated by the coefficient matrix in 
RREF 

/I 1 1 \ 
A= 1 -1 . 

yo 1 2 / 



7 



-50 50 100 -100 -50 50 100 

Figure 2.2: Line soliton solutions of KPII: (a) the (3,2)-soliton solution generated by the coefficient matrix A in 
Eq. (12. 8b with (^i,...,^5) = (— 1,0,|,|,|) at t = —32; (b) the (3,3)-soliton solution generated by the coefficient 
mattix A in Eq. ^ with , . . . jA;^) = (-1, -^,0, 5, 1, j) at f = 20. 

Again, we first determine the asymptotic line solitons as y ^ 00. In this case, the pivot columns of A are 

labeled by {£1,^2,^3} = {1,4,5}. So, the asymptotic line solitons as y —s- 00 are given by the index pairs 

[1,71], [4,72] and [5,73] for some ji, 72, 73. Starting with the first pivot, ei = 1, we take j = 2,3,... and 

/I \ 

check the rank of the sub-matrix Xu in each case. When / = 2 we have Xu = 1 -1 , so rankfXi?) = 

Vo 1 2 / 

3 > A'^ — 1. Hence, according to Propositio ri2.4f ii). the index pair [1,2] does not correspond to an asymptotic 

/I \ 

line soliton as y ^ 00. We then take 7 = 3 and consider the sub-matrix Z13 = ^0 -1 j . Since rank(Xi3) = 2 

and rank(Xi3|l) = rank(Xi3|3) = rank(Xi3[l,3) = 3, the rank conditions in Propositio ri2.4f ii) are satisfied. 
Therefore the index pair [1 , 3] corresponds to an asymptotic line soliton as 3^ ^ 00. Moreover, by considering 
7 =4,5,6 one can easily check that the rank conditions are no longer satisfied. Thus [1,3] is the unique 
asymptotic line soliton associated with the pivot index = 1 as y — > 00. Proceeding in a similar way, we find 

that for the pivot column e? = 4, the rank conditions are only satisfied when 7 = 5, since X45 = ( -1 ) , is 

Vo 2 / 

of rank 2, and rank(Z45|4) = rank(X45|5) = rank(X45|4, 5) = 3. Therefore, the index pair [4,5] corresponds 
to an asymptotic line soliton as y ^ 00. Finally, we find that the third asymptotic line soliton as y ^ 0° is 
given by the index pair [5,6]. 

Next, we proceed to determine the asymptotic line solitons as y — >^ —00. The non-pivot columns of A 
are labeled by the indices g\ = 2 g2 = ?> and ^3 = 6. For g\ = 2, the only possible value of / < 7 is / = 1. 
In this case F12 = 0, so rank(7i2) = and rank(Fi2ll) = rank(Fi2|2) = rank(7i2| 1,2) = 1. Thus, the pair 
[1,2] identifies an asymptotic line soliton as y ^ —00. For ^2 = 3 we consider / = 2, 1. When / = 2, the 
rank conditions in Proposition I2.41 ii) are satisfied, leading to the asymptotic line soliton [2,3] as y — > —00. 
Similarly, it is easy to verify that for ^3 = 6 the index pair [4, 6] uniquely identifies the asymptotic line soliton 
as y ^ —00. Summarizing, there are = 3 asymptotic line solitons as y — > 00 identified by the index pairs 
[1,3], [4,5] and [5,6], and there are = 3 asymptotic line solitons as y ^ —00 identified by the index pairs 
[1,2], [2,3] and [4,6]. A snapshot of this (3,3)-soliton solution at f = —20 is shown in Fig. 12.2b . 

So far we have discussed the properties of the generic (A^-jA/^^) -soli ton solutions of KPII, and have 
shown how to characterize the asymptotic line solitons as y — > ±00 from the corresponding T-function. In 
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the next section, we investigate an important subclass of the (A^-jA/^^) -soli ton solutions called the elastic 
A^-soliton solutions. 

3 Elastic A/^-soliton solutions 

We begin this section by introducing the notion of equivalence classes of soliton solutions and their T- 
functions, both of which will play important roles in this section. 

Definition 3.1 (Equivalence class) Let denote the set of all exponential phase combinations whose coef- 
ficients are non-zero in the z-function ofEq. (12.41 ). Two tau-functions are said to be in the same equivalence 
class if they contain the same set @ (up to an overall exponential phase factor). All {N-,N+)-soliton so- 
lutions of KPII generated by an equivalence class of T-functions form an equivalence class of solutions. 

It is clear from the above definition that the T-functions in a given equivalence class can be viewed as 
positive-definite sums of the same exponential phase combinations but with different sets of coefficients. 
Such T-functions are parametrized by the same set of phase parameters ki,... ,kM, but the constants 0„,o in 
the phase 6^ are different. Moreover, the irreducible coefficient matrices associated with the T-functions 
have exactly the same sets of vanishing and non-vanishing minors, but the magnitudes of the non-vanishing 
minors are different for different matrices. Then, an important consequence of Proposition I2.2l is that for all 
(A^-jA^^^) -soli ton solutions in the same equivalence class, the corresponding asymptotic line solitons arise 
from the same i —>■ j single-phase transition, and are therefore labeled by the same index pair [i,j]. Propo- 
sition |23] then implies that the coefficient matrices associated with the T-functions in the same equivalence 
class have identical sets of pivot and non-pivot indices which identify respectively, the asymptotic line soU- 
tons as J — > oo and as j ^ — oo. Thus, solutions in the same equivalence class can differ only in the position of 
each asymptotic line solitons and in the location of each interaction vertex. As a result, any (A^_,A'^+) -soli ton 
solution of KPII can be transformed into any other solution in the same equivalence class by spatio-temporal 
translations of the individual asymptotic line solitons. 

The KPII equation (11.11) is invariant under the symmetry {x,y,t) — s- {—x,—y,—t). Thus, if u{x,y,t) 
is an (M — A'^jA'^) -soliton solution of KPII, then u{—x,—y,—t) is also an (A^_,A'^+) -soliton solution whose 
incoming and outgoing line solitons are reversed so that A'^- = A^ and A'^-i- = M — N. In general, u{x,y,t) 
and u{—x,—y,—t) belong to two different equivalence classes of solutions, and so do their generating T- 
functions. However, the function T{—x,—y,—t) generating the solution u{—x,—y,—t) is not by itself a 
T-function according to Eq. (12.41 ). Using Eq. (12.41 ). ■T{—x,—y,—t) can be expressed as T{—x,—y,—t) = 
i'(x,y,t) , where the function 

z'{x,y,t)= A{mi,...,mN) exp[d{h,...,lM-N)] Yi (^"v"^"'.) 

l<mi<-<mjv<M l<.v<r<A' 

is a positive definite sum of exponential phase combinations labeled by the set of indices {l\, . . . Jm-n}, 
which is the complement of {m\, . . . ,mN} in {1,2,...,M}. Moreover, since T'{x,y,t) only differs from 
x{—x, —y, —t) by an overall exponential phase factor, it should be clear from Eq. (12.11) that they both generate 
the same solution u{—x, —y, —t). The correspondence between equivalence classes of solutions and their t- 
functions related via the symmetry {x,y,t) {—x, —y, —t) leads to the following notion of duality. 
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Definition 3.2 (Duality) Two equivalence classes of x -functions are said to be dual if they are parametrized 
by the same set of phase parameters {^1,^2, • • • j'^m}- but correspond to complementary sets of exponential 
phase combinations and &. That is, 

e{mi,...,mN) €@ <^ 0(/i,...,/m-w) G0', 

where the sets {mi, . . . ,mi^} and {h, - ■ ■ ,Im-n} form a disjoint partition of the integers {1,2, .. . ,M}. Sim- 
ilarly, two equivalence classes of {N-,N^)-soliton solutions are dual if they are generated by dual equiva- 
lence classes of z-functions. 

In particular, if a given (M — A'^,A'^)-soliton solution u{x,y,t) belongs to a certain equivalence class, then the 
corresponding (A/^jM — A'^)-soliton solution u{—x, —y, —t) belongs to the dual equivalence class. 

An interesting subclass of (A^-jA/^^) -soli ton solutions are the elastic A/^-soliton solutions of KPII as men- 
tioned in section [U These can be defined as follows. 

Definition 3.3 A {N-,Nj^)-soliton solution is called elastic if it belongs to an equivalence class which is its 
own dual. 

Clearly, in this case we have A'^+ = = A'^ and M = 2N. Moreover, the amplitudes and directions of the 
A'^ incoming line solitons coincide with those of the A'^ outgoing line solitons. Thus, an elastic A'^-soliton 
solution is generated by a "self-dual" T-function which is a positive definite sum over a set of exponential 
phase combinations such that the following condition holds: 

e(mi,...,mw) G0 ^ 0(/i,...,/;v) G0, V{mi, . . . ,mw} U {/j , . . . = {1,2, . . . ,2A^} . 

We outline below the main properties of the elastic A'^-sohton solutions. Additional details regarding these 
solutions can be found in Refs. BlfTTl. 

Proposition 3.4 The elastic N-soliton solutions of KPII are characterized as follows. 

(i) Each elastic N-soliton solution u{x,y,t) and the dual solution u{—x, —y, —t) belong to the same equiv- 
alence class. 

( ii) The T-function corresponding an elastic N-soliton solution has M = IN distinct phase parameters and 
anN X 2N, irreducible, rank N coefficient matrix A whose N xN minors satisfy the duality conditions: 

A{mi,...,mN)=0 ^ A{Ii,...,In) = 0, (3.1) 

where the indices {mi , mf^} and {/i , . . . , lp)}form a disjoint partition of integers {1,2,..., 2A'^}. 

C Hi) Each elastic N-soliton solution exactly N asymptotic line solitons as y ^ ±00 identified by the same 
index pairs with e„ < g„, n = 1, . . . ,A'^. The indices e\,e2^. ■ ■ ^e^^ and gi,g2, ■ ■ ■ ,gN label re- 

spectively, the pivot and non-pivot columns of the coefficient matrix A. Hence, they form a disjoint 
partition of integers {1,2,..., 2A^}. 

(iv) The amplitude and direction of the asymptotic line soliton [e„,g„] are the same as y ^ ±00, and 
are given in terms of the phase parameters as an = kg^ — kg^ and Cn = kg^ + k^,^. 
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The set Sn := {(a„,c„)|fl'„ > 0}^^^ C M.^^ of all admissible A/^-tuples of amplitude-direction pair associ- 
ated with elastic A/^-soliton solution will be called the soliton parameter space. An element {a„,c„}^^j G 
of soliton parameters is admissible, if it yields a set of 2N distinct phase parameters Kp^ := {kf}^^y where 
= (c« i fl'n)/2 and <k^ ■ By sorting the elements of Kj^ in increasing order k\ < .. . < kzN, one obtains 
the ordered set = {^i ,^2; • • • ,^2^} of phase parameters associated with the T-function. The positions of 
the phase parameters k^ of the line soliton can be labeled uniquely within the ordered set K by an ordered 
pair of indices [inj,,] such that /„ < j„. That is, k^ = ki^, k^ = kj^^. Since [/„,7„] also identify two distinct 
columns of the coefficient matrix A, it follows that /„ = g„ which labels a pivot column, and j„ = gn which 
labels a non-pivot column of A, for elastic A/^-soliton solutions. We refer to this identification between the 
sets Kn and K as phase pairing which defines a map Sn — > S, where S is the set of all possible choices of 
N distinct integer pairs {[in-,jn]}^=i from {1,2, .. . ,2N}. This map identifies each set of soliton parameters 
{ian,Cn)}n=i & Sn to & sct {[e,i,gn]}^^^i G 5 of (pivot, non-pivot) index pairs. Note however that distinct 
elements of Sn can in fact lead to the same pairing if the elements of the corresponding sets Kn are ordered 
in exactly identical fashion. Thus, phase pairing induces a partition of the soliton parameter space Sn into 
several disjoint sectors. Each sector is distinguished by a single element {[^,i,g,j]}^^i G 5 of distinct integer 
pairs which labels the A'^ asymptotic line soli tons corresponding to any set {('3m>c„)}^^j of soliton param- 
eters chosen from that sector in Sn- Therefore, the total number of such disjoint sectors is given by the 
number of elements of the set 5 namely, l^l = {2N — 1) ! !. Furthermore, given any set of soliton pai^ameters 
{(<3„,c„)}^^j from one of these sectors in Sn, it is possible to construct a coefficient matrix A in RREF 
satisfying Conditions 12.11 and 12.31 and whose pivot and non-pivot columns are labeled by the corresponding 
set of index pairs {[^n,gn]}^=i G S obtained via phase pairing. The matrix A is constructed by using the 
rank conditions of Proposition I2.4f ii) and the duality condition Eq. (13.11 ) in Proposition I3.4f ii). However, it 
is not unique but is determined up to some free parameters. Therefore, the matrix A obtained in this way 
together with the set K of phase parameters, produce an equivalence class of T-f unctions from (12.41 ). The 
directions and amplitudes of the asymptotic line soUtons in the corresponding equivalence class of elastic 
A/^-soliton solutions coincide with the set of soliton parameters {(fl'n,c„)}^^j G Sn that was originally chosen. 
We illustrate these facts in the following example where we explicitly construct an elastic 3-soliton solution 
only from its soliton parameters. 

Example: We start with the set {(ai ,ci), (fl2,C2), (a3,C3)} = {(5/2, —5/2), (5/2, — 1/2), (7/4,5/4)} G S3. 
From this we construct the set of unordered phase parameters ^"3 := {^^ = (c„ ± £?„) /2},"j^j whose elements 
are kf = 0, A:f = -5/2, k:^ = \,k:^ = -3/2, )t+ = 3/2, k:^ = -1/4. Note that the elements of K3 are all 
distinct. Sorting these elements of ^^3 in increasing order, we obtain the ordered set of phase parameters 
K = {—5/2, —3/2, — 1/4, 0, 1 , 3/2}. Comparing the elements of the sets Ki, and K, we obtain the following 
phase pairing: k^ = k\, k^ = k4, k2 = k2, k2 = k(,, kj = kj,, k^ = k^. Hence, the asymptotic line solitons 
are labeled by the index pairs [1,4], [2,6] and [3,5], where {e\,e2,es} = {1,2,3} are the pivot indices and 
{§1,^25^3} = {4,6,5} are the non-pivot indices of the corresponding coefficient matrix A. Note also that 
the pivot indices are sorted but the non-pivot indices are unsorted. Next, we outline the construction of 
the coefficient matrix A in RREF satisfying Conditions 12. 1 1 and 12.31 and whose pivot and non-pivot column 
indices are specified above. The construction proceeds in several steps. 

Step 1. We start with a 3 x 6 matrix in RREF with pivot and non-pivot columns labeled respectively, by the 
indices {1,2,3} and {4,6,5}: 

(1 Ml Vl Wl \ 
1 -U2 -V2 -W2 , 
1 M3 V3 W3 / 
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where {«/, v/, w/}^^j are nonnegative numbers. The negative signs in the second row arise due to Condi- 
tion [ZTJc) which demands that if any given 3x3 minor of A is non-zero, then it must be positive. 

Step 2. In order to obtain further information about A, we apply the rank conditions in Proposition I2.4r ii) 
to the sub-matrices and F,y associated with each line soliton [/, j]. For example, starting with the line 

soliton [1,4] as J — > — oo and by considering the sub-matrix F14 = (10), we find that rank(Fi4) = 2. Then, 

rank(Fi4|4) must be 3, which means that the minor A(2, 3,4) = mi 7^ 0. Now suppose vi = 0, then the non- 
negativity of the minors in Condition 12. U c) implies that A(2,4,5) = —u\Vt, > and A(3,4,5) = —u\V2 > 0, 
whose only solution is V2 = V3 = since mi 7^ 0. But then the 5* column of A contains no nonzero elements, 
violating the irreducibility of A in Condition 12.31 Thus, vi 7^ 0, and similar arguments lead to wi 7^ 0. Then 
from A(/, a, j8) > where / G {2,3} and a,j3 € {4,5,6}, we can deduce that if ui = 0, then v/ = w/ = 
for I = 2,3. Consequently, the only nonzero element in each of the 2"'* and 3'^'' row of A would be the pivot 
entry, and again this would violate Condition [231 Hence, we must also have M2 7^ and M3 7^ 0. As our goal 
is to obtain only one representative matrix A associated to the equivalence class of T-function, we simplify 
subsequent calculations by choosing a particular normalization such that nonzero elements mi = vi = wi = 1. 
Then we have 

/I 1 1 1 \ 

A = 1 -U2 —V2 -W2 ■ 

\0 1 M3 V3 W3 / 

10 1 

Step 3. Next, we consider the line soliton [3,5] as y 00 and the associated sub-matrix X35 = (01 -W2 ) . 



From the condition rank(X35) < N — I = 2, we have det(X35) = A(l,2,6) =0, which implies that 1V3 = 0. 
Moreover, since the minor A(l,2,6) = 0, it follows from the duality condition Eq. (13.11 ) that A(3,4,5) = 
U2 — V2 = 0. Hence V2 = M2 7^ 0. Applying the duality condition again to the minor consisting of the 
pivot columns, we obtain A(l,2,3) = 1 7^ ^ A(4,5,6) 7^ 0. In particular, this means that the 5* and 

6"^ columns of A are linearly independent, and that the sub-matrix Xu = ("V2 -W2) associated with the 
[1,4] line soliton as j — > 00, has rank 2. Then it follows from the rank conditions in Proposition I2.4f ii) 
that rank(Xi4|l) =3 ^ A(l,5,6) = V3 W2 7^ 0. Thus, we have V3 7^ and W2 7^ 0. Finally, imposing the 
non-negativity condition on the remaining minors, we obtain the following form of the coefficient matrix A 

/I 1 1 1 \ 

A = 1 —U2 —U2 —W2 

\0 1 M3 V3 / 

where the remaining free parameters satisfy < W2 < M2 and < V3 < M3. Thus, starting only from 
the soliton parameters, we have constructed a 4-parameter family of coefficient matrices corresponding 
to an equivalence class of elastic 3-soliton solutions whose asymptotic line sohtons are labeled by the 
index pairs [1,4], [2,5], [3,6]. An elastic 3-soliton solution generated by the above coefficient matrix A 
with {u2,W2,ut,,vt,) = (1,2/3,2/3,3/5) and K= {—5/2, —3/2, —1/4,0, 1,3/2} (as above), is shown in 
Fig. [He). 

By further investigating the combinatorial properties of the coefficient matrix A, it is possible to obtain 
additional information regarding the classification scheme for the elastic A^-soliton configuration space. 
These results are presented below. 

Proposition 3.5 Each elastic N -soliton configuration is described by a set {[en-,8n\\^=i of distinct integer 
pairs with e„ <gn,n= I,. .. ,N. The indices e„ label the pivot columns and the indices gn label the non-pivot 
columns of the irreducible coefficient matrix A in RREF. In addition, the following results hold: 
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(i) Without any loss of generality, the pivot indices can be ordered as \= e\ < e2 < ■ ■ ■ < efj < 2N. Each 
pivot index e„n= \ . ,N satisfies the inequality n < e„ <2n — I. Moreover, the number of possible 

(2N)\ 

ways of choosing the pivots are given by the Catalan number Cn = jy-j^jy-_|_ ^-^ j - 

(ii) For a fixed ordered set of pivots, the number of possible choices of a (unordered) set 

N 

{gi ) • • • jSn} of non-pivot indices such that e„ < V« = 1, . . . ,A^, is given by ]^(2?i — 

n=l 

( Hi) The total number of ways of choosing N distinct pairs in the set S is given by [IN — 1 ) ! !. 



Note that the requirement that the set of integer pairs {[en,gn\]^^=\ be distinct was already stated (without 
proof) in Ref. ifTTI . and some of the above-hsted consequences were also obtained there. 

We illustrate the results in Proposition 13.51 by presenting the classification scheme for the elastic 3- 
soliton solution space. This is achieved by enumerating all possible arrangements of the pivot positions in 
the irreducible coefficient matrix A in RREF. In this case, N = 3 and A is a 3 x 6 matrix with 3 pivots whose 
possible (column) positions are determined by Proposition 13. 5f i) as follows: = 1, 2 < ^2 < 3, 3 < ^3 < 5. 
Thus, the total number of pivot configurations is given by C3 = 6!/(3!4!) = 5. Thus, this classification 
scheme gives rise to 5 subclasses of elastic 3-soliton solutions. The number of inequivalent types of solutions 
in each subclass is determined by all possible {[eM)^n]}^=i pairings for a given choice of the pivot positions 
{ei,e2,eT,}. These are obtained from Proposition I3.5f ii). and are itemized below. 

3 

(i) Pivot positions: {^1,^21^3} = {li2,3}. Total number of distinct pairings = ]^(2?i — e„) = 3! = 6. 

«=i 

List of inequivalent elastic 3-soliton solutions: 

{[1,4], [2,5], [3,6]}, {[1,5], [2,4], [3,6]}, {[1,6], [2,4], [3,5]}, 
{[1,4], [2,6], [3,5]}, {[1,5], [2,6], [3,5]}, {[1,6], [2,5], [3,4]}. 



(ii) Pivot positions: {^1,^2,^3} = {1,2,4}. Total number of distinct pairings = 4. List of inequivalent 
elastic 3-soliton solutions: 

{[1,3], [2,5], [4,6]}, {[1,3], [2,6], [4,5]}, {[1,5], [2,3], [4,6]}, {[1,6], [2,3], [4,5]}. 



(iii) Pivot positions: {^1,^2,^3} = {1,2,5}. Total number of distinct pairings = 2. List of inequivalent 
elastic 3-soliton solutions: {[1,3], [2,4], [5,6]}, {[1,4], [2,3], [5,6]}. 

(iv) Pivot positions: {£1,^2,^3} = {1,3,4}. Total number of distinct pairings = 2. List of inequivalent 
elastic 3-soliton solutions: {[1,2], [3,5], [4,6]}, {[1,2], [3,6], [4,5]}. 

(v) Pivot positions: {^1,^2, ^3} = {1,3,5}. Total number of distinct pairings = 1. Elastic 3-soliton solu- 
tion: {[1,2], [3,4], [5,6]}. 

Thus, the total number of inequivalent 3-soliton solutions is6 + 4 + 2 + 2 + l = 15 = 5!!as given by Propo- 
sition [33tiii). Fig. 13. II shows a sample from the fifteen inequivalent cases. Fig. l3.1f a) shows the previously 
known ordinary 3-soliton solution (cf. section [T]), while the remaining solutions are new, and they exhibit 
resonant interactions. 
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(a) [1,2], [3,4J, L5,6J (b) [1,3], [2,5], [4,6] (c) [1,4], [2,6], [3,5] (d) [1,5], [2,4], [3,6] 



Figure 3.1: Examples of elastic 3-soliton solutions. 

4 Conclusion 

In this article, a family of multi-soliton solutions of the KPII equation have been studied. These solutions are 
generated by T-functions which are expressed as positive definite linear combinations of exponential phases 
that are linear in the variables x,y,t. It is remarkable that such a simple form of the T-function generates 
multi-soliton configurations which exhibit a rich variety of time dependent spatial structures including res- 
onant interactions and web patterns. The asymptotic analysis of the faw-function in the xy-plane reveals that 
the solution decays exponentially except along certain directions which are characterized by the transition 
between two dominant exponential phase combinations which have all but one phases in common. All such 
("non-decaying") directions for any given solution can be explicitly identified by analyzing the x M ma- 
trix coefficient A associated with the T-function. In particular, as j ^ oo there are directions which can be 
identified with the pivot columns of A; while as y — > — oo, there are M — N directions which can be identified 
with the non-pivot columns of A. 

When M = 2N, the general line soliton solutions contain the special subclass of the elastic //-soliton 
solutions. Each elastic A^-soUton solution has a set of directions in the xy-plane as y — > oo, and an identical 
set of A^ directions as y —>■ — oo along which the solution does not decay. Moreover, for any given A^, the 
solution space of elastic A'^-solitons can be decomposed into {2N — 1)!! distinct regions, each region cor- 
responding to inequivalent types of solutions. It is interesting to note that the previously known ordinary 
A^-soUtons form only one of these types. Thus, the space of elastic A^-soUton solutions of KPII appears to be 
much richer than previously thought. 

It is significant that solutions exhibiting siinilar features of soliton resonance and web structure have 
also been obtained in several other (2+l)-dimensional integrable systems, besides KPII. These solutions 
were also derived by direct algebraic methods similar to the approach taken here. Therefore, it is reasonable 
to expect that the results developed in this work for KPII will also be useful to characterize soliton solutions 
in a variety of other (2+l)-dimensional integrable systems. 
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